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Abstract 

We propose the measurement of the ratio Rcss(Qr m ) = a ^ T> ® T ^ to study 
the effects of the multiple soft gluon radiation, predicted by QCD, on the 
transverse momentum (Qt) distribution of the weak gauge bosons and Z° 
produced at the Tevatron. We compare the prediction of the extended Collins- 
Soper-Sterman resummation formalism with the next-to-leading and next-to- 
next-to-leading order calculations. We show that both the rich dynamics of 
the QCD multiple soft gluon radiation and the non-perturbative sector of 
QCD can be tested by measuring Rcss- 
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I. INTRODUCTION 



With 100 ptT 1 luminosity at the Tevatron, about 2 x 10 6 W ± and 6 x 10 5 Z° bosons 
are produced, and the data sample will increase by a factor of 20 in the Run 2 era. In 
view of this large event rate, a careful study of the transverse momentum distributions of 
vector bosons can provide a stringent test of the rich dynamics of the multiple soft gluon 
emission predicted by Quantum Chromodynamics (QCD). The increasing precision of the 
experimental data demands a high precision theoretical calculation of the distributions of the 
W ± and Z° bosons, which takes the effects of multiple gluon radiation into account. In this 
work, within the extended Collins-Soper-Sterman formalism [|I]-|3[], we illustrate the effect 
of the multiple gluon radiation on the vector boson transverse momentum distribution. 

To test the dynamics of the multiple soft gluon radiation, in this work we propose the 
measurement of the ratio Rcss(Qr m ) = CT ^J r > ^ ; T - for the W ± and Z° bosons produced at 
the Tevatron. We show that for vector boson transverse momenta less than about 30 GeV, 
the difference between the resummed and the fixed order predictions (either at the as or 
a s order) can be observed by measurement. This suggests that in this kinematic region, 
the effects of the multiple soft gluon radiation are important, and hence, this case provides 
a ideal opportunity to test this aspect of the QCD dynamics. For Qt less than about 10 
GeV, the Qt distribution is largely determined by the non-perturbative part of QCD. At the 
Tevatron, for and Z° production, this non-perturbative physics, when parametrized by 
Eq. (U), is dominated by the parameter g2, which was shown to be related to properties of the 
QCD vacuum 0]. Therefore, precisely measuring the Qt distributions in the low Qt region, 
e.g. from Z° events can advance our knowledge of the non-perturbative QCD physics. 

The next section summarizes the relevant formulae of our extension of the Collins-Soper- 
Sterman resummation formalism, which describes the production and decay of vector bosons 
at hadron colliders. In Section fTJ, we present the differences in the next-to-leading order 
(NLO, e.g. 0(azs)), next-to-next-to-leading order (NNLO, e.g. 0(a s )), and the resummed 
predictions for Rcss, and show that the differences due to the soft gluon effects are mea- 
surable. We illustrate that, using the experimental data, the improvement of the non- 
perturbative sector of the resummation formalism is also possible. Finally, we draw our 
conclusions, based upon these theoretical results, about the importance of the multiple soft 
gluon radiation in and Z° production at the Tevatron. 
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II. SUMMARY OF THE RESUMMATION FORMALISM 



The fully differential cross section of the hadronic production and decay of a vector 
boson, within the extended Collins-Soper-Sterman resummation formalism, is characterized 
as follows. The kinematics of the vector boson V (real or virtual) can be expressed in the 
terms of its mass Q, rapidity y, transverse momentum Qt, and azimuthal angle </y, measured 
in the laboratory frame (the center-of-mass frame of hadrons h\ and h 2 ). The kinematics 
of the lepton t x is described by 9 and 0, the polar and the azimuthal angles defined in the 
Collins-Soper frame 0, which is a special rest frame of the \^-boson ||. (A more detailed 
discussion of the kinematics can be found in Ref. ||.) The resummed cross section is given 
by the following formula in Ref. 0: 

fda(h 1 h 2 ^V(->£ 1 I 2 )X)\ 1 Q 2 



dQ 2 dy dQ 2 T d<f) V d cos 6d<f)J res 96tt 2 S {Q 2 - My) 2 + Q 4 T 2 V /M 2 

-i- / d 2 be i ^ rb J2Wf k (K,Q,x 1 ,x 2 ,9,cj ) ,C 1 ,C 2 ,C 3 )Wf k p (b,Q,x 1 ,x^ 

+ Y(Q T ,Q,x 1 ,x 2 ,9,<f>,C i )}. (1) 

In the above equation the parton momentum fractions are defined as x\ = e y Q/\fS and 
x 2 = e~ v Q/y/S, where \/~S is the center-of-mass (CM) energy of the hadrons hi and h 2 . 
For V = or Z°, we adopt the LEP line-shape prescription of the resonance behavior. 
The renormalization group invariant quantity Wj k (b), which sums to all orders in a g all the 
singular terms that behave as asQj> 2 hi 2 " 1 " 1 (Qt/Q 2 ) (I < rn < n) for Qt — > 0, is 

W jk (b, Q, xi, x 2 , 9, <f>, C h C 2 , C 3 ) = exp {-5(6, Q, C h C 2 )} | V jk \ 2 
x { [(C ja ® f a/hl ) (an) (C- kb ® f b/h2 ) (x 2 ) + (C~ ka ® f a/hl ) (xx) (c jb ® f b/h2 ) (x 2 ) 
x(9l + 9 2 R )(f L + f R )(l + cos 2 9) 

+ 

x 



C ja ® /a/hx) (Xl) (C kb ® f b/h2 ) (X 2 ) - (C ka <g> fa/hi) (»l) (p jb ® f b /h 2 ) (x 2 ) 

(9 2 L-9 2 n)(f 2 L-f 2 n)(^os9)}, (2) 



where <S> denotes the convolution 

(C ja ® / a/hl ) (xi) = J ^ C ja 6, ^ = C x , C 2 ^) /a/ft! (fi, A* = ^ 

and the coefficients are given by 

J Cabibbo — Kobayashi — Maskawa matrix elements for V = W ± 
Vjk = { 6 jk for V = Z° 
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The qq'V and l\l 2 V vertices are defined as ry^c^l — 75) + gn(l + 75)] and ^[/^(l — 75) + 
fn(l + 75)], respectively. For example, for V = W + , q — u, q' = d, l\ = v e , and I 2 = e + , 
the couplings are g\ = j\ — GpM^/\/2 and g 2 R = }r = 0, where Gf is the Fermi constant. 
The Sudakov exponent S{b, Q, Ci, C 2 ) in Eq. (g) is defined as 



S(b,Q,C 1 ,C 2 )= f 



r,2 



Cl/b 2 /I 2 



A (a s (fi), d) In f^M + ^ (a s (/x), d, C 2 ) 



The explicit forms of the A, i? and C functions and the renormalization constants Ci (z=l,2,3) 
are summarized in Refs. 

In Eq. ([]]) the magnitude of the impact parameter b is integrated from to 00. However, 
in the region where b 3> 1/Aqcd, the Sudakov exponent S(b, Q, Ci, C 2 ) diverges as the 
result of the Landau pole of the QCD coupling as{[i) at fi = Aqcd, and the perturbative 
calculation is no longer reliable. In this region of the impact parameter space (i.e. large 
b), a prescription for parametrizing the non-perturbative physics in the low Qt region is 
necessary. Following the idea of Collins and Soper 0, the renormalization group invariant 
quantity Wjj.{b) is written as 

W fk (b) = W fk {K)W% p (b) . 

Here Wjj~(b*) is the perturbative part of Wjj.{b) and can be reliably calculated by perturbative 
expansions, while W^ p {b) is the non-perturbative part of Wjj,{b) that cannot be calculated 
by perturbative methods and has to be determined from experimental data. To test this 
assumption, one should verify that there exists a universal functional form for this non- 
perturbative function W^ p (b). This is similar to the general expectation that there exists a 
universal set of parton distribution functions (PDF's) that can be used in any perturbative 
QCD calculation to compare it with experimental data. In the perturbative part of Wjj,(b), 

h 



V 1 + 



max I 



and the non-perturbative function was parametrized by (cf. Ref. JTJ) 



w fk P ( b i Q> Qo, X!, x 2 ) = exp 



-F 1 (b) In - F j/hl (x l7 b) - F i/ha (x 2 , b) 



0, 



(3) 



where Fx, Fj/ hl and h^/ h2 have to be first determined using some sets of data, and later can 
be used to predict the other sets of data to test the dynamics of multiple gluon radiation 
predicted by this model of the QCD theory calculation. As noted in Ref. Jl|], Fi does not 
depend on the momentum fraction variables x\ or x 2 , while Fj/^ and Fj.i h2 in general depend 
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on those kinematic variables.[| The \b.(Q 2 /Qq) dependence associated with the F± function 
was predicted by the renormalization group analysis Furthermore, Fi was shown to 
be universal, and its leading behavior (~ b 2 ) can be described by renormalon physics Q. 
Various sets of fits to these non-perturbative functions can be found in Refs. || and ||. 

In our numerical calculations, we use the Ladinsky-Yuan parametrization of the non- 
perturbative function (cf. Ref. [Bf): 



Wf k p {b,Q,Q ,x 1 ,x 2 ) =exp ,,, ... 



9ib 2 - g 2 b 2 ln ( ) - gig 3 b In (100xix 2 ) 



(4) 



where g x = O.llt&os GeV 2 , # 2 = 0.58±&£ GeV 2 , # 3 = -1.5tjft GeV" 1 , and Q = 1-6 GeV. 
(The value 6 ma:r = 0.5 GeV -1 was used in determining the above g^s.) These values were 
fit for CTEQ2M PDF with the canonical choice of the renormalization constants, i.e. C\ = 
C*3 = 2e _TB (7£ is the Euler constant) and C 2 = 1. In principle, for a calculation using a 
different set of PDF, these non-perturbative parameters should be refit using a data set that 
should include the recent high statistics Z° data from the Tevatron. 

In Eq. (P, Wji sums over the soft gluon contributions that grow as a^Q^ 2 hi 2 ™ -1 (Qt/Q 2 ) 
(1 < m < n) to all orders in as- Contributions less singular than those included in Wj% 
should be calculated order- by-order in as and included in the Y term, introduced in Eq. (|l|). 
This would in principle extend the applicability of the CSS resummation formalism to all 
values of Qt-Q The Y term, which is defined as the difference between the fixed order 
perturbative contribution and those obtained by expanding the perturbative part of Wj% to 
the same order, is given by 



Y(Q T , Q, Xl , x 2 , 9, 0, C 4 ) = f 1 *j± f f 2 - f; 



a s {C 4 Q) 
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xfaMti, C A Q) R ( ^(Qt, Q,^,^,6,4>) c*Q), (5) 

where the functions contain contributions less sing ular than agQ^ 2 In 2 ™" 1 (Q 2 T /Q 2 ) 
(1 < m < n) as Qt — ► 0. Their explicit expressions are summarized in Refs. Pfl. 



1 Here, and and throughout this work, the flavor dependence of the non-perturbative functions is 
ignored, as it is postulated in Ref. Q. 

2 It is shown in Ref. [fj] that since the A, B, C, and Y functions are only calculated to some finite 
order in as, the CSS resummed formula as described above will cease to be adequate when the 
value of Qt is in the vicinity of Q. Hence, in practice, one has to switch from the resummed 
prediction to the fixed order perturbative calculation as Qt > Q- 
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III. THE RATIO Rcss 



In this work we propose to measure the ratio Rcss(Qr m ) = ^(Qt > Qt" 1 ) 1 1 a Total to 
distinguish the predictions of the resummed, NLO and NNLO calculations. In Fig. [1] we 
show the distributions of Rcss, which is defined by 

p ^mim _ °(Qt > Qf n ) 1 da(hh 2 ^V) 
McssWt > = = / «VT 



OTotal CfTotal JQ^ in " ' dQ T 

where Q™ 8 * is the largest Qt allowed by the phase space. In the NLO calculation, a(Qr > 
Qt 1u ) grows without bound near Q™ m = 0, as the result of the singular behavior 1/Q T in the 
matrix element. The NLO curve runs well under the resummed one in the 2 GeV < Q™ in < 
30 GeV region, and the Qt distributions from the NLO and the resummed calculations have 
different shapes even in the region where Qt is of the order of 15 GeV. 

With a large number of fully reconstructed Z° events at the Tevatron, one should be able 
to use the data to clearly discriminate these two theoretical calculations. The experimental 
uncertainty in the total cross sections of the W ± and Z° productions, based on 19.7 pb _1 
CDF data, is in the ballpark of 5% |IJ. Fig. |l] shows that, in the 10 GeV < Q^' m < 30 



GeV region, even with this experimental precision we should see deviations between the 
experiment and the NLO predictions, in which the effects of the multiple gluon radiation 
are not included. In view of this result it is not surprising that the DO analysis of the as 



measurement [11] based on the measurement of a(W + 1 jet)/a(W + jet) does not support 
the NLO calculation. We expect that if this measurement were performed demanding the 
transverse momentum of the jet to be larger than about 50 GeV, at which scale the resummed 
and NLO distributions cross (cf. Ref. |§), the NLO calculation would adequately describe 
the data. 

To show that for Qt below 30 GeV the QCD multiple soft gluon radiation is important 
to explain the DO data [[U]] , we also include in Fig. [I] the prediction for the Rcss distribution 
at the order of a s . As shown in the figure, the a s curve is closer to the resummed curve, 
which proves that for this range of Qt the soft gluon effect included in the a s calculation is 
important for predicting the vector boson Qt distribution. In other words, in this range of 
Qt, it is more likely that soft gluons accompany the boson than just a single hard jet 
associated with the vector boson production. For large Qt, it becomes more likely to have 
hard jet(s) produced with the vector boson. 

The measurement of Rcss can a ls° provide information about the non-perturbative 
physics associated with the initial state hadrons. As shown in Ref. ||, the effect of the 
non-perturbative physics on the Qt distributions of the W ± and Z° bosons produced at the 
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FIG. 1. The ratio Rcss as a function of Q™ m for W + bosons. The fixed order [0(ots) short 
dashed, 0{a 2 s ) dashed] curves are ill-defined in the low Qt region. The corresponding distributions 
for the Z° boson are indistinguishable from those for the in this plot. 

Tevatron is important for Qt less than about 10 GeV. This is evident if one observes that 
different parametrizations of the non-perturbative functions do not change the Qt distri- 
bution for Q T > 10 GeV, but can dramatically change its shape for Q T < 10 GeV. Since 
for W ± and Z° production, the Ib.(Q 2 /Qq) term is large, the non-perturbative function, as 
denned in Eq. (|3|), is dominated by the F\(b) term which is expected to be universal for all 
Drell-Yan type processes and which is related to the physics of the renormalon ||]. Hence, 
the measurement of Rcss can be used to probe this part of non-perturbative physics for 
Qt < 10 GeV, in addition to probing the dynamics of multiple soft gluon radiation for 10 
GeV < Qt < 40 GeV. It is therefore important to measure Rcss at the Tevatron. With a 
large sample of Z° data at the Run 2, it will be possible to determine the dominant non- 
perturbative function which can then be used to determine the W ± boson Qt distribution to 
improve the accuracy of the My/ and the charged lepton rapidity asymmetry measurements. 

IV. CONCLUSIONS 

In conclusion, the measurement of Rcss provides an accurate test of the dynamics of the 
multiple soft gluon radiation predicted by QCD. With high enough luminosity, the NLO, 
NNLO, and resummed theoretical predictions can be distinguished by the and Z° data 
at the Tevatron. The comparison of the NNLO and resummed predictions shows that the 
soft gluon effect is important in the Qt < 30 GeV region. Additionally, in the Qt < 
10 GeV region, the Rcss measurement, using the recent Z° data at the Tevatron, can 
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provide valuable information on the non-perturbative sector of the resummation formalism. 
Therefore, a careful measurement of the Q T distribution of the vector bosons W ± and Z° 
in the Q T < Q/2 region at the Tevatron can further our knowledge of the perturbative 
dynamics and the non-perturbative domain of QCD. 
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